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The geometrically frustrated spin-1/2 Ising-Heisenberg model on triangulated Husimi lattices is 
exactly solved by combining the generalized star-triangle transformation with the method of exact 
recursion relations. The ground-state and finite-temperature phase diagrams are rigorously calcu¬ 
lated along with both sublattice magnetizations of the Ising and Heisenberg spins. It is evidenced 
that the Ising-Heisenberg model on triangulated Husimi lattices with two or three inter-connected 
triangles-in-triangles units displays in a highly frustrated region a quantum disorder irrespective of 
temperature, whereas the same model on triangulated Husimi lattices with a greater connectivity of 
triangles-in-triangles units exhibits at low enough temperatures an outstanding quantum order due 
to the order-by-disorder mechanism. The quantum reduction of both sublattice magnetizations in 
the peculiar quantum ordered state gradually diminishes with increasing the coordination number 
of underlying Husimi lattice. 
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I. INTRODUCTION 

Over the last few decades, two-dimensional (2D) frus¬ 
trated Heisenberg models have become subject of in¬ 
tense research studies owing to an exceptional diver¬ 
sity in a low-temperature magnetic behavior they dis¬ 
play [H-Q. However, the rigorous treatment of 2D frus¬ 
trated Heisenberg models is mostly unattainable due to 
extraordinary difficulties, which bear a close relation to 
non-commutative (quantum) character of spin operators. 
Many semi-classical 2D Ising spin systems with compet¬ 
ing interactions may also exhibit intriguing cooperative 
phenomena originating from a geometric spin frustration, 
but their experimental realizations are regrettably scarce 
[^, @ at the expense of exact solvability of the respec¬ 
tive models [ 3 , 0 ] • Generally, the spin frustration is in¬ 
dispensable ground for emergence of several remarkable 
cooperative phenomena such as the order-by-disorder, 
which refers to a selection of unusual ordered state of 
some highly frustrated classical or quantum spin model 
driven either by thermal or quantum fluctuations [Ml. 
Another striking phenomenon closely related to the ge¬ 
ometric spin frustration is a reentrance, which occurs if 
an increase in temperature causes a phase transition from 
a disordered state emerging at lower temperatures to a 
spontaneously ordered state emerging at higher temper¬ 
atures Q. The reentrant phase transitions have been 
experimentally detected in various magnetic systems in¬ 
cluding spin glasses [l3l| or even non-magnetic systems 
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such as liquid mixtures M- 

Exactly solved models belong to the most attracting 
issues of statistical mechanics, because they provide re¬ 
sults not affected by any approximation [15j. It is well 
known that the exact solution of the Ising model on the 
Bethe lattice can alternatively be viewed as the rather re¬ 
liable (Bethe-Peierls) approximation of this model on a 
regular planar lattice with the same coordination number 
pU [r^ . However, the Ising model on the Bethe lattice 
cannot capture consequences of a geometric spin frus¬ 
tration inherent to the antiferromagnetic Ising model on 
close-packed planar lattices, because the Bethe lattices 
do not contain any (odd-numbered) cycles. This defi¬ 
ciency can be avoided by considering the Ising model on 
another type of recursive lattice built up from the inter¬ 
connected polygons, whereas the interior part deep inside 
of such recursively built tree is usually referred to as the 
Husimi lattice [isl - l^ . It is worthy to notice that the 
Ising and Ising-like models defined on the Husimi lat¬ 
tices have played an important role in a description of 
several interesting physical systems and phenomena like 
binary alloys, rare gases, lipid b ilay ers, spin liquids, spin 
glasses or polymers (see Ref. and references cited 
therein). 

Recently, the spin-1/2 Ising-Heisenberg model on tri¬ 
angulated (triangles-in-triangles) planar lattices have 
attracted a great deal of attention [29l - l33| . because 
it closely resembles a magnetic structure of real 
polymeric coordination compounds Cu 9 X 2 (cpa) 6 -nH 20 
(X=E. Cl,Br) with the geometry of triangulated kagome 
lattice |26l - l^ . Besides a clear experimental motivation, 
the spin-1/2 Ising-Heisenberg model on triangulated pla¬ 
nar lattices has afforded a valuable playground for an 
investigation of the role of local quantum fluctuations 
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within a novel class of exactly solvable classical-quantum 
spin systems. It has been convincingly evidenced (3lj | 
that the spin-1/2 Ising-Heisenberg model with a greater 
connectivity of triangles-in-triangles units in the under¬ 
lying triangular lattice gives rise to stronger local quan¬ 
tum fluctuations, which preferably support a remarkable 
quantum order at low enough temperatures instead of a 
quantum disorder through the order-by-disorder effect. 
Contrary to this, the spin-1/2 Ising-Heisenberg model on 
triangulated kagome lattice cannot exhibit this intriguing 
feature [2^ [s^ . Bearing this in mind, it appears worth¬ 
while to examine the spin-1/2 Ising-Heisenberg model on 
related triangulated Husimi lattices, in which the con¬ 
nectivity of triangles-in-triangles units can be varied sys¬ 
tematically. 

The outline of this paper is as follows. The spin-1/2 
Ising-Heisenberg model on triangulated Husimi lattices 
will be introduced in Sec. |IT] together with a few basic 
steps of its exact treatment. The most interesting re¬ 
sults for the ground-state and finite-temperature phase 
diagrams will be discussed along with typical thermal 
variations of spontaneous magnetization in Sec. IHII The 
paper ends up with a brief summary of the most impor¬ 
tant findings in Sec. IIVI 


II. MODEL AND METHOD 


natively can be rewritten as a sum over the cluster Hamil¬ 
tonians 

Nq/3 

H ^ Hfc, (2) 

fe=i 

whereas each cluster Hamiltonian TLk involves all the 

interaction terms of three Heisenberg spins from the 

kth triangle-in-triangle unit (i.e. the six-spin cluster 
schematically illustrated at the bottom of Fig. [T]) 

3 

3 

= (3) 

i=l 

It is quite evident that the standard commutation re¬ 
lations = 0 are satisfied for different cluster 

Hamiltonians i ^ j and hence, the partition function 
of the spin-1/2 Ising-Heisenberg model on triangulated 
Husimi lattices can be partially factorized into the fol¬ 
lowing product 

Nq/3 Nq/3 

-Z^IHM = En Trfeexp(-/3?/fc) = En Zk, (4) 

{o-i} fc=l {o-i} fc=l 


Let us consider the spin-1/2 Ising-Heisenberg model 
on triangulated Husimi lattices, which are schematically 
illustrated on the left-hand-side of Fig. [1] for two se¬ 
lected values of the coordination number q = S and 4. 
The overall magnetic structure of triangulated Husimi 
lattices form smaller triangles of three quantum Heisen¬ 
berg spins (Heisenberg trimers), which are embedded into 
each larger triangle of the pure Husimi lattice involving 
in its nodal lattice sites the classical Ising spins. The to¬ 
tal Hamiltonian of the spin-1/2 Ising-Heisenberg model 
on triangulated Husimi lattice denoting a deep interior 
such recursively built tree can be defined as follows 


where /3 = l/(fcBT), /cb is the Boltzmann’s constant, T 
is the absolute temperature, the symbol X]{cri} denotes a 
summation over spin states of all the Ising spins and the 
symbol Tr^ = Trsj,jTrs^ 2 Trs ,^3 stands for a trace over 
spin degrees of freedom of the kth Heisenberg trimer. 
The mathematical structure of Eq. (|4]) implies that one 
may perform a trace over spin degrees of freedom of dif¬ 
ferent Heisenberg trimers independently of each other 
in order to obtain an explicit form of the Boltzmann’s 
weight Zk , which will solely depend on three nodal Ising 
spins tTfci, (Tfc 2 , and Cfea attached to the kth Heisenberg 
trimer. Afterwards, one may take advantage of the gen¬ 
eralized star-triangle transformation 


Nq 2Nq 

n = -JnJ2+ stsf) + -JiJ2 

{k,l) {k,j) 

( 1 ) 

Here, and {a = x, y, z) label spatial components of 
the usual spin-I/2 operator, the parameter Jh accounts 
for the XXZ Heisenberg interaction between the nearest- 
neighbor Heisenberg spins, A is a spatial anisotropy in 
this interaction term, the parameter Ji stands for the 
Ising interaction between the nearest-neighbor Heisen¬ 
berg and Ising spins, respectively, N denotes the total 
number of the Ising spins and q is the coordination num¬ 
ber that determines how many triangles-in-triangles units 
meet at each nodal site of the underlying Husimi lat¬ 
tice. The total Hamiltonian o of the spin-I/2 Ising- 
Heisenberg model on triangulated Husimi lattices alter¬ 


Zk (o-fei, cr^ 2 : = Trfe exp(-/3?ife) 

= Aexp[/3R{alia^2 + <^k2<^k3 + '^k3<^ki)]^ ( 5 ) 


which replaces the effective Boltzmann’s weight Zk with 
the equivalent expression depending on three nodal Ising 
spins only. The star-triangle mapping transformation 
generally represents a set of eight algebraic equations, 
which can be obtained from Eq. (O by substituting all 
available spin conhgurations of three Ising spins involved 
therein. In fact, one gets just two independent equations; 
the first equation for two uniform configurations with the 
three equally aligned Ising spins 


Zi = 2 exp ( -13 Jh I cosh ( -/3Ji 


-I- 2 exp 


-i/3JH(l-4A) 


cosh 
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FIG. 1: The spin-1/2 Ising-Heisenberg model on the triangulated Husimi lattices (figures on the left) and its rigorous mapping 
equivalence with the effective spin-1/2 Ising model on the pure Husimi lattices (figures on the right) established by means of the 
generalized star-triangle transformation (figure at the bottom). The triangulated Husimi lattices, which include two different 
lattice sites occupied by the Ising (empty circles) and Heisenberg (full circles) spins, are drawn for two particular values of the 
coordination number q = 3 and 4. 


-I- 4 exp 
= A exp 


--/3Jh(1 + 2A) 


cosh 1 -pji ) , 




( 6 ) 


and the second one for six nonuniform configurations 
with one Ising spin pointing in opposite with respect to 
the other two 


1 


Z 2 = 2 exp I -/3Jh cosh -/3Ji 


+ 2 exp 
-I- 2 exp 
-I- 2 exp 


--/3Jh(1 + 2A) 

-i/3JH(l-A) 

-i/3JH(l-A) 


1 


cosh ( -/3Ji 
cosh ( 


cosh ( 


Aexp ( -^I3R ) , 


(7) 


which involves the parameters defined as follows 
Q"*" = + 2 (JrA)^. 


V 2 


( 8 ) 


It is noteworthy that the algebraic equations m and © 
unambiguously determine yet unspecified mapping pa¬ 


rameters A and R 

A={Z^Zl)^'\ = (9) 

which enter into the generalized star-triangle transforma¬ 
tion ([ 5 ]). On assumption that the mapping parameters 
A and R satisfy Eqs. dSl)-®, the star-triangle trans¬ 
formation m holds quite generally (regardless of spin 
states of three nodal Ising spins) and it can be substi¬ 
tuted into the relation (U) with the aim to get a rigorous 
mapping correspondence between the partition function 
ZmM of the spin-1/2 Ising-Heisenberg model on the tri¬ 
angulated Husimi lattices and, respectively, the partition 
function Zim of the effective spin-1/2 Ising model on a 
pure Husimi lattice 

2ihm(/3, Jh, Ji, A) = A^'^/^ZimW. R), (10) 

which is unambiguously given by the Hamiltonian with 
the effective nearest-neighbor Ising interaction R: 

Nq 

Him =—7? ^ erf cr|. (H) 

(hi) 

It should be noted that the effective spin-1/2 Ising 
model on pure Husimi lattices can be exactly solved by 
means of the method of exact recursion relations, which 
consequently enables to extract exact results for the spin- 
1/2 Ising-Heisenberg model on the triangulated Husimi 
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lattices from the relevant mapping correspondence be¬ 
tween both models. The rigorous treatment of the spin- 
1/2 Ising model on triangular Husimi lattices using the 
recursion method has been reported on previously in sev¬ 
eral foregoing works [isl - l^ . so it is sufficient to recall 
the most important outcomes of this calculation. It can 
be easily proved that the effective nearest-neighbor cou¬ 
pling dni) of the corresponding spin-1/2 Ising model on 
pure Husimi lattices is always ferromagnetic (i.e. R > 0), 
which substantially simplifies the iteration procedure as 
three inter-connected recursion relations can be replaced 
with a single recursion relation. Assuming the ferro¬ 
magnetic pair interaction R > 0 between the nearest- 
neighbor spins, the spontaneous magnetization of the 
spin-1/2 Ising model on triangular Husimi lattices with 
the coordination number q can be expressed as 

mjM = (o-fci)iM = 2 ’ (12) 

where the symbol (■ ■ piM denotes a canonical ensem¬ 
ble average performed within the effective Ising model 
defined by the Hamiltonian dill) and x represents the 
stable fixed point of the recurrence relation (i.e. x = 
llcUyi—^oo ^n) 

1 + -f exp(/3i?)x^'(ri^ 

eMm) + 2xl-_\+xl^ir 

Similarly, the expectation values of other important 
quantities such as the triplet correlation function between 
three Ising spins residing sites of an elementary trian¬ 
gular plaquette of the pure Husimi lattice can also be 
calculated from the stable fixed point of the recurrence 
relation (fT3l) according to the formula 


Husimi lattices from the corresponding rigorous results of 
the effective spin-1/2 Ising model on triangular Husimi 
lattices. It can be straightforwardly proved by adopt¬ 
ing the exact mapping theorems developed by Barry et 
al. [srl - li^ that the sublattice magnetization of the Ising 
spins within the spin-1/2 Ising-Heisenberg model on tri¬ 
angulated Husimi lattices directly equals to the sponta¬ 
neous magnetization of the effective spin-1/2 Ising model 
on pure Husimi lattices 


= (<5’fci)lHM = (o’^i)lM = TOIm(/3A). 


(17) 


Here, the symbol (• • pmM denotes a canonical ensemble 
average performed within the spin-1/2 Ising-Heisenberg 
model on the triangulated Husimi lattice defined through 
the Hamiltonian o and the single-site magnetization 
of the effective spin-1/2 Ising model on the triangular 
Husimi lattices can be calculated from Eqs. m and 
dT51) after solving the latter recursion relation iteratively. 
Contrary to this, the sublattice magnetization of the 
Heisenberg spins in the spin-1/2 Ising-Heisenberg model 
on triangulated Husimi lattices can be exactly calculated 
from the generalized form of Callen-Suzuki spin identity 

EMI 


mn 



/ Trfc;ggiexp(-/3?7fc) \ 
\ Trfeexp(-/3-Hfc) / 


According to Eq. (TTSl) . the sublattice magnetization of 
the Heisenberg spins mn can be related to the formerly 
derived single-site magnetization of the Ising spins mi 
given by Eqs. (fT^ - (fT^ and the triplet correlation func¬ 
tion tiM given by Eqs. dH-dlSl) 


miM 

Qi Q 2 

2 tiM 

Qi 

-3^ 

2 

Zi Z 2 

3 

Zi 

^2 


(19) 


Hm — (0'fcl0’fe2^fc3)lM ~ 8 ( \ + X'i ) ’ 


which involves another two parameters y and z defined 
through the stable fixed point of recurrence relation (fT^ 

exp(/3i?) — 2x‘‘~^ + 

^ exp(/3i?) -|- 2x^~^ + ’ 

1 — 2x'^~^ + exp(/?i?)a;^'^“^ 
exp(/3i?) -I- 2x‘i~'^ x‘^i~‘^ 

It has been rigorously proved by Jurcisinova and Jurcisin 
[ 2 ^ that the spin-1/2 Ising model on triangular Husimi 
lattices with the ferromagnetic nearest-neighbor interac¬ 
tion R > 0 exhibits a continuous (second-order) phase 
transition from a ferromagnetic phase to a paramagnetic 
phase at the following (inverse) critical temperature 


PcR 


R 

k-eTc 


In 


2q + l 
2g- 3 


(16) 


Two new functions Qi and Q 2 entering the exact formula 
m for the sublattice magnetization of the Heisenberg 
spins are defined as follows 


Qi = 3 exp I -/3 Jh I cosh I -/3Ji 


exp -i/3JH(l-4A) sinh Q/3 Ji 


2 exp --/3 Jh (1-I-2A) 


1 


sinh 


Q 2 = 3 exp -/3Jh sinh -/3Ji 


exp 

exp 

exp 


-i/3JH(l + 2A) 

-i/3JH(l-A) 

-i/3JH(l-A) 


sinh 1 -/3Ji 


cosh 


cosh ( ^ 1^*5 


( 20 ) 


With all this in mind, let us extract exact results for 
the spin-1/2 Ising-Heisenberg model on the triangulated 


Last but not least, the critical points of the spin-1/2 
Ising-Heisenberg model on triangulated Husimi lattices 


























5 



FIG. 2: (Color online) The ground-state phase-diagram of the 
spin-1/2 Ising-Heisenberg model on the triangulated Husimi 
lattices in the A — Jh/Ji plane. 


can be readily obtained from a comparison of the ef¬ 
fective coupling with its critical value m, which 
subsequently provides the following explicit form for the 
critical condition 

{2q-3)Z1 = {2q + l)Zl (21) 

where the superscript c means that the inverse criti¬ 
cal temperature /3c enters into the expressions Zi and 
Z 2 defined by Eqs. ® and © instead of the usual 
inverse temperature /3. The critical condition m 
can be rather straightforwardly used for constructing 
finite-temperature phase diagrams of the spin-1/2 Ising- 
Heisenberg model on triangulated Husimi lattices with 
arbitrary value of the coordination number q. 


III. RESULTS AND DISCUSSION 

Let us proceed to a discussion of the most interesting 
results for the spin-1/2 Ising-Heisenberg model on tri¬ 
angulated Husimi lattices. In the following, we will be 
mainly interested in addressing the question of whether 
the coordination number of the underlying Husimi lattice 
may essentially influence the overall magnetic behavior of 
the studied model or not. In general, the change of the 
ferromagnetic Ising interaction (Ji >0) to the antiferro¬ 
magnetic one (Ji < 0) does not cause any fundamental 
difference in a respective magnetic behavior and hence, 
our further attention will be merely restricted to a spe¬ 
cial case with the ferromagnetic Ising coupling Ji > 0 
henceforth serving as the energy unit. By contrast, the 
antiferromagnetic Heisenberg interaction Jh < 0 is re¬ 
sponsible for a mutual interplay of local quantum fluc¬ 
tuations with geometric spin frustration totally absent 
in the model with the ferromagnetic Heisenberg interac¬ 
tion Jh > 0 , which makes an investigation of the highly 
frustrated region Jh/Ji ^ 0 particularly intriguing. 

Let us examine first a ground state of the spin-1/2 
Ising-Heisenberg model on triangulated Husimi lattices. 


It is worth mentioning that the overall ground-state 
eigenvector is given by a tensor product over the lowest- 
energy eigenstates of the cluster Hamiltonians ([3]), which 
means that the ground-state phase diagram and individ¬ 
ual ground states of the spin-1/2 Ising-Heisenberg model 
on triangulated Husimi lattices are completely identical 
with that reported on previously for the analogous model 
on triangulated planar lattices composed of the same 
structural unit |2iM^ . In this regard, let us merely 
quote individual ground-state eigenvectors of the spin- 
1/2 Ising-Heisenberg model on triangulated Husimi lat¬ 
tices without their detailed description (the interested 
reader is referred to Ref. for further details). Alto¬ 
gether, the ground-state phase diagram depicted in Fig. 
[3] in the A — Jn/Ji plane includes two unique sponta¬ 
neously ordered ground states, namely, the classical fer¬ 
romagnetic (CF) ground state 

N Nq/3 

|CF) = nit>< n imj.,.;,,.;,. (22) 

2=1 k —1 


the quantum ferromagnetic (QF) ground state 

N Nq/3 

I QF)= n Itlr n (Itn)+int)+.5^3(23) 


and one highly degenerate ground-state manifold 
spanned over the following 10 eigenstates of the cluster 
Hamiltonians ([3]) 




^<2 

■'^k2 




;)5(im)-iin))s^3 


'^k2 


’^k2 


’^k2 




’3k2 




'3k2 


^3^k2 


'3k2 


^3k2 



(24) 


It should be noted that the highly degenerate ground- 
state manifold (1331) basically changes its character in the 
limiting Ising case A = 0, which does not exhibit any 
interesting feature due to a complete lack of local quan¬ 
tum fluctuations and this special case will be therefore 
left out from our subsequent ground-state considerations. 
Although the spin-1/2 Ising-Heisenberg model on most 
triangulated planar lattices [2^-[3l| prefers a quantum 
disorder (QD) in a highly frustrated parameter region 
pertinent to the ground-state manifold (1331) . it has been 
convincingly evidenced in Ref. 3l| that the model alter¬ 
natively may display a spectacular spontaneous quantum 
order (QO) due to the order-by-disorder effect. These 
results would suggest that the geometry of underlying 
magnetic lattice can basically influence an essence of the 
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ground-state manifold (l24l) . which is usually disordered 
but under certain conditions may become spontaneously 
ordered. 

Let us find a rigorous criterion that will definitely 
resolve the puzzling issue concerning with the sponta¬ 
neously ordered or disordered nature of the ground-state 
manifold (0^ of the spin-1/2 Ising-Heisenberg model on 
triangulated Husimi lattices. It is quite evident from Eqs. 
(El and p9l) that a nonzero spontaneous magnetization 
of the effective spin-1/2 Ising model on the pure Husimi 
lattices represents sought condition ensuring nonzero 
spontaneous magnetizations of the Ising and Heisenberg 
spins. The spin-I/2 Ising model on the triangular Husimi 
lattices has nonzero spontaneous magnetization if and 
only if the effective nearest-neighbor coupling /3i? exceeds 
the critical value given by Eq. (HU), which monotoni- 
cally decreases with rising coordination number of the 
underlying Husimi lattice (e.g. /3cI? = ln(5) for q = 2, 
PcR = ln(7/3) for q = 3, /3ci? = ln(9/5) for q = A, 
etc.). In the highly frustrated region Jh/Ji < —2-1 (2 +A) 
relevant to the ground-state manifold the effective 
coupling acquires the zero-temperature asymptotic value 
limT_>o /3i? = In 2 that is conclusive whether or not a 
spontaneous long-range order appears. A comparison of 
the zero-temperature asymptotic limit lim 7 ’_>.o /3i? = In 2 
with the critical value (HSl) furnishes a rigorous proof 
that the ground-state manifold (IMl) remains disordered 
just for the spin-1/2 Ising-Heisenberg model on the tri¬ 
angulated Husimi lattices with the coordination numbers 
q = 2 and 3, while the same model on the triangulated 
Husimi lattices with greater values of the coordination 
number g > 4 shows a peculiar spontaneous long-range 
order of quantum nature. 

To verify our theoretical prediction of the quantum 
order-by-disorder phenomenon, let us explore in detail 
finite-temperature phase diagrams of the spin-1/2 Ising- 
Heisenberg model on the triangulated Husimi lattices, 
which are plotted in Fig. [3] for three different values of 
the exchange anisotropy A and a few values of the coordi¬ 
nation number q. First, it can be seen from Fig. [Sj/a) that 
the spin-1/2 Ising model on triangulated Husimi lattices 
retrieved in the limit A = 0 shows a familiar monotonous 
dependence of the critical temperature on the interaction 
ratio Jh/<3i, which becomes zero whenever the antiferro¬ 
magnetic interaction Jh/</i < —1 is strong enough to 
force the investigated model to the disordered ground 
state. The highest value of the critical temperature is 
accordingly obtained in the asymptotic limit Jh/^i oo 


lim 

Jt-f 


3 Jl 



I + 2^2(9-1) 
2q-3 



(25) 


On the contrary, the spin-I/2 Ising-Heisenberg model 
on the triangulated Husimi lattices exhibits more diverse 
finite-temperature phase diagrams due to a mutual in¬ 
terplay between the geometric spin frustration and lo¬ 
cal quantum fluctuations, which are only present if the 
anisotropy parameter A > 0. Fig. [2Kb) illustrates vari¬ 





FIG. 3: (Color online) The critical temperature as a function 
of the ratio between the Heisenberg and Ising interaction for 
several values of the coordination number q and three different 
values of the exchange anisotropy: (a) A = 0; (b) A = 1; and 
(c) A = 2. 


ations of the critical temperature with the interaction 
ratio Jh/Ti that are quite typical for the model with 
the isotropic Heisenberg coupling A = 1 (the displayed 
case) as well as the XXZ Heisenberg coupling of the 
easy-axis type 0 < A < 1. Both aforementioned cases 
indeed show qualitatively the same dependences except 
a small quantitative difference including the asymptotic 
limit Jh/Ji —f oo, at which the critical temperature still 
reaches the asymptotic value (ESI) in the case of easy-axis 
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anisotropy 0 < A < 1 but it goes to a slightly lower 
asymptotic value in the isotropic case A = 1. How¬ 
ever, the most important hnding relates to the critical 
frontiers in a highly frustrated region. In accordance 
with our ground-state argumentation, the critical tem¬ 
perature tends to zero at the ground-state boundary 
(Jh/Ji = —2/3 for A = 1) between the classical ferro¬ 
magnetic phase and the highly degenerate ground- 
state manifold (l24l) just for the spin-1/2 Ising-Heisenberg 
model on two triangulated Husimi lattices with the coor¬ 
dination numbers q = 2 and 3. Note furthermore that the 
model dehned on the triangulated Husimi lattice with the 
coordination number g = 3 displays a subtle reentrant re¬ 
gion when the interaction ratio is selected slightly below 
the relevant ground-state phase boundary Jh/^/i ^ —2/3. 
On the other hand, the critical temperature of the spin- 
1/2 Ising-Heisenberg model on triangulated Husimi lat¬ 
tices with greater coordination numbers (? > 4 remains 
nonzero in the whole frustrated region (Jh/Ji < —2/3 
for A = 1), which provides an independent confirmation 
of the spectacular quantum order arising out from the 
highly degenerate ground-state manifold (1^ due to the 
quantum order-by-disorder effect. The greater the co¬ 
ordination number of the triangulated Husimi lattice is, 
the higher is the critical temperature of the spontaneous 
quantum long-range order. 

The last hnite-temperature phase diagram shown in 
Fig. [Sic) illustrates the critical temperature of the spin- 
1/2 Ising-Heisenberg model on the triangulated Husimi 
lattices when assuming the XXZ Heisenberg coupling of 
the easy-plane type A > 1. Two appropriate comments 
are in order. First, it is quite apparent from Fig. |3Kb) 
and Fig. EIc) that the critical temperature qualitatively 
retains character of its dependence on the interaction ra¬ 
tio Jh/Ji in the highly frustrated region for any A > 0, 
i.e. regardless of whether the Heisenberg interaction has 
the easy-axis anisotropy, the easy-plane anisotropy or 
is completely isotropic. Second, the critical tempera¬ 
ture unexpectedly shows a considerable decrease upon 
strengthening of the ferromagnetic Heisenberg coupling 
in a nonfrustrated regime Jh/Ji > 0. The intriguing sup¬ 
pression of the critical temperature can be attributed to 
a phase transition from the classical ferromagnetic phase 
(1^^ to the quantum ferromagnetic phase (051) . which oc¬ 
curs at Jh/Ji = 1 for the selected value of the exchange 
anisotropy A = 2. Owing to this fact, the critical tem¬ 
perature approaches the asymptotic limit 


lim 


ksTc 

Ji 


1 f, \2q-l + 2^2{q-l) 

3 f' [-2jV3- 



(26) 


which is for the quantum ferromagnetic phase exactly 
three times lower than the one (l25ll of the classical fer¬ 
romagnetic phase on account of a quantum reduction of 
the spontaneous magnetization of the Heisenberg spins to 
one-third of its saturation value. Hence, it follows that 
the quantum entanglement and associated quantum fluc¬ 
tuations are responsible for a gradual suppression of the 





FIG. 4; (Color online) Thermal variations of the sublattice 
magnetizations of the Ising spins (solid lines) and the Heisen¬ 
berg spins (broken lines) for several values of the interaction 
ratio Jh/Ji and a few values of the coordination number and 
the exchange anisotropy: (a) 5 = 3, A = 1; (b) q = 3, A — 2; 
(c) g = 4, A= 1. 


critical temperature of the quantum ferromagnetic phase 
(1231) in response to a strengthening of the ferromagnetic 
Heisenberg interaction of easy-plane type. 

Last but not least, let us discuss typical temperature 
changes of spontaneous sublattice magnetizations of the 
Ising and Heisenberg spins, respectively. It can be de¬ 
duced from Fig. 0] that both sublattice magnetizations 
usually follow the same general trends though the sub¬ 
lattice magnetization of the Heisenberg spins generally 








































TABLE I: Zero-temperature asymptotic limits of the sublat¬ 
tice magnetization of the Ising and Heisenberg spins as a func¬ 
tion of the coordination number q of the triangulated Husimi 
lattices in the highly frustrated region Ju/Ji < ~2/(2 -|- A). 


4 

mi 

m-H 

g < 3 

0 

0 

g = 4 

0.3543 

0.1394 

g = 5 

0.4537 

0.1630 

g = 6 

0.4808 

0.1659 




q ^ OO 

1/2 

1/6 


falls down more rapidly with an increase in temperature 
than the sublattice magnetization of the Ising spins. It 
can be seen from Fig. Hl^a) that the critical temperature 
of the spin-1/2 Ising-Heisenberg model on the triangu¬ 
lated Husimi lattice with the coordination number q = 3 
gradually diminishes as the antiferromagnetic Heisenberg 
interaction strengthens owing to reinforcement of the ge¬ 
ometric spin frustration. The most notable tempera¬ 
ture dependence of both sublattice magnetizations can 
be detected for a sufhciently but not too strong anti¬ 
ferromagnetic Heisenberg interaction (Jh/Ti < —2/3 for 
A = 1), which drives the investigated model to the dis¬ 
ordered ground state (1241) and simultaneously preserv¬ 
ing it in a close vicinity of the phase boundary with 
the classical ferromagnetic phase (|2^ . Under these cir¬ 
cumstances, the sublattice magnetizations of the spin- 
1/2 Ising-Heisenberg model on the triangulated Husimi 
lattice with the coordination number g = 3 display 
reentrance with two consecutive critical temperatures in 
concordance with the finite-temperature phase diagram 
shown in Fig. |3Jb) (see in Fig. HJa) the dependence for 
Jh/Ji = -0.68). 

To confirm an existence of the quantum ferromagnetic 
phase (1^ we have depicted in Fig. SKb) temperature 
variations of both sublattice magnetizations of the spin- 
1/2 Ising-Heisenberg model on the triangulated Husimi 
lattices with the ferromagnetic Heisenberg interaction of 
the easy-plane type. While the sublattice magnetization 
of the Ising spins remains fully saturated irrespective of 
a relative strength of the ferromagnetic Heisenberg in¬ 
teraction, the sublattice magnetization of the Heisenberg 
spins is fully saturated just if Jh/Ji < 1 when consider¬ 
ing the special case with the anisotropy parameter A = 2. 
Assuming the reverse condition Jh/Ji > 1 and A = 2, 
the sublattice magnetization of the Heisenberg spins un¬ 
dergoes a quantum reduction to one-third of its satura¬ 
tion value due to a symmetric quantum superposition of 
three up-up-down spin states of the Heisenberg trimers. 
This finding evidently corroborates an existence of the 
quantum ferromagnetic ground state (1^51) in overwhelm¬ 
ing parameter space with the ferromagnetic Heisenberg 
coupling of the easy-plane type. 

Of course, the most challenging is validation of the 


predicted spectacular quantum order in a highly frus¬ 
trated region of the spin-1/2 Ising-Heisenberg model on 
the triangulated Husimi lattices with sufficiently high co¬ 
ordination numbers g > 4. To illustrate the case, Fig. 
SKc) shows typical temperature dependences of both sub¬ 
lattice magnetizations of the spin-1/2 Ising-Heisenberg 
model on the triangulated Husimi lattice with the coor¬ 
dination number g = 4 when considering the isotropic 
Heisenberg interaction A = 1. It can be readily veri¬ 
fied that the sublattice magnetizations of Heisenberg and 
Ising spins are subject to a quantum reduction caused 
by local quantum fluctuations. This assertion is con¬ 
sistent with the observed zero-temperature limits of the 
sublattice magnetizations, which acquire in the highly 
frustrated region {Ja/Ji < —2/3 for A = 1) nontriv¬ 
ial asymptotic values. Tab. U summarizes the zero- 
temperature limits of the sublattice magnetization of 
Ising and Heisenberg spins for a few selected values of 
the coordination number g, which were obtained by solv¬ 
ing a polynomial equation of degree 2g — 1 derived from 
the exact recursion relation (fT51) . It can be understood 
from Tab. Uthat the quantum reduction of both sublat¬ 
tice magnetizations gradually shrinks as the coordination 
number g of the triangulated Husimi lattice increases. 
The quantum reduction of the spontaneous magnetiza¬ 
tion of the Ising spins completely disappears in the limit 
g —^ oo in contrast to the spontaneous magnetization of 
the Heisenberg spins, which tends towards one-third of its 
saturation value in the limit g —> oo. Finally, the notable 
temperature-induced increase in both sublattice magne¬ 
tizations can be observed in Fig. SKc) for Jh/Ji ^ —2/3, 
which bears a close relation to a thermal excitation from 
the ground-state manifold (1^ to the classical ferromag¬ 
netic phase (1221) accompanied with an increase in the rel¬ 
evant magnetizations (see the curves for Jh/Ji = —0.68 
and —0.75). 


IV. CONCLUSIONS 

The present work deals with the exact solution of 
the spin-1/2 Ising-Heisenberg model on the triangulated 
Husimi lattices, which has been obtained by the gen¬ 
eralized star-triangle mapping transformation and the 
method of exact recursion relations. In particular, our 
attention has been focused on a rigorous determination of 
the ground-state and finite-temperature phase diagrams, 
which were subsequently corroborated by temperature 
dependences of the sublattice magnetizations of Ising and 
Heisenberg spins serving as the relevant order parame¬ 
ters. It has been convincingly evidenced that the spin- 
1 /2 Ising-Heisenberg model defined on two closely related 
triangulated Husimi lattices may display fundamentally 
different magnetic behavior even if both lattices may dif¬ 
fer from each other just in a connectivity of the same 
structural triangles-in-triangles unit. 

We have furnished rigorous proof that the quantum 
disorder develops in a highly frustrated region of the spin- 
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1/2 Ising-Heisenberg model on the triangulated Husimi 
lattices with the coordination number q = 2 and 3, while 
the identical model on the triangulated Husimi lattices 
with greater coordination numbers <? > 4 exhibits the 
spectacular quantum order due to the order-by-disorder 
effect. In addition, it has been demonstrated that the 
spontaneous magnetizations of the Ising and Heisenberg 
spins undergo a quantum reduction in the remarkable 
quantum ordered state, which gradually shrinks with in¬ 
creasing a connectivity of the triangles-in-triangles units. 
It should be also mentioned that reentrant phase tran¬ 


sitions can only be found in the model defined on the 
triangulated Husimi lattice with the coordination num¬ 
ber (7 = 3. To conclude, the spin-1/2 Ising-Heisenberg 
model on the triangulated Husimi lattices represents a 
novel exactly solved classical-quantum spin model, which 
has proved its usefulness in elucidating a sought connec¬ 
tion between the quantum order-by-disorder effect and 
the relative strength of local quantum fluctuations var¬ 
ied systematically through the coordination number of 
the underlying Husimi lattice. 
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